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ABSTRACT 

This study was undertaken to investigate the effects 
of varying both the secondary Mach nuitiber and angle of 
inclination relative to the primary stream on the flow 
field generated by the interaction of a secondary jet with 
a supersonic mainstream. The experimental portions of 
this investigation were conducted at a primary Mach number 
of 2.80 in the Naval Postgraduate School Supersonic Wind 
Tunnel. Data are presented and compared with various 
theories. This presentation includes correlation of the 
penetration height of the secondary flow, bow-shock shape 
by non-dimensionalization with respect to various param- 
eters and correlation of interaction side force amplifi- 



cation factor. 
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I. 



INTRODUCTION 



For a period of time spanning at least the last decade, 
considerable interest has been displayed in the interaction 
of secondary gaseous flows injected into a primary super- 
sonic stream. Several theoretical models exist which 
reflect the physical trends expected for varying sets 
of injection parameters (such as injection pressure, 
temperature, mass flow rate, Mach number and angle of incli 
nation to the primary flow, et cetera) . The body of experi 
mental investigation has, for the most part, concerned 
itself with sonic injection normal to the primary stream. 

The objective of this work is to study the phenomena 
associated with changes of both injection Mach number and 
inclination angle in order to compare results with the 
analytical models presently available. The degree to which 
correlation of experimental and theoretical results are 
obtained and the empirical relationships observed will be 
discussed. 

This study was performed at the Naval Postgraduate 
School Supersonic Blow-Down Wind Tunnel operating at a 
primary Mach number of 2.80. A flat plate model incor- 
porating an interchangeable nozzle set and instrumented 
for the collection of plate surface pressure distribution 
was constructed and used in this investigation. Commer- 
cially available dry nitrogen, helium and argon were 
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injected into the main air flow at various stagnation 
pressures and mass flow rates through nozzles of varying 
Mach number and angle of inclination to the primary 
stream. 
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II. EXPERIMENTAL APPARATUS 



A. THE FLAT PLATE MODEL 

The pylon mounted flat plate model shown in Figure 1 
(made of steel and flash plated with chromium) was mounted 
in the test section of the Naval Postgraduate School 
Supersonic Wind Tunnel. The tubing which leads to the 
forty-seven pressure taps arranged on the upper surface 
of the flat plate passes through the model pylon and the 
tunnel floor. The pressure tap locations are depicted 
in Figure 2. Access for mating the nozzles with the 
injectant gas supply lines is also made available through 
the model pylon and tunnel floor. A representative set of 
drawings of three interchangeable injection nozzles is 
provided in Figure 3 showing, from left to right, a normal 
sonic nozzle, a normal supersonic nozzle and an inclined 
supersonic nozzle. All nozzles have a throat diameter of 
0.1000 ± 0.001 and convergence and divergence half-angles 
of 10°. Nozzles were constructed for injection at Mach 
numbers of 1.00, 1.73 and 2.13. For each Mach number 
nozzles were constructed at nominal inclination angles e 
of 90°, 80°, 70°, 60° and 50°. Epsilon is taken to be the 
angle between the nozzle axis and the direction of the 
primary stream (or the surface of the flat plate) in a 
sense such that for injection not normal to the primary 
flow an upstream component of the secondary flow exists. 
The nozzles are inserted from above the model. Alignment 



15 



and a positive in-place lock for securing the nozzles in 
position is provided by a small key bolt passing through 
the pylon from below 

B. THE WIND TUNNEL 

This investigation was conducted in the Naval Post- 
graduate School Blow-Down Wind Tunnel. The facility has a 
test section 4 inches by 4 inches in cross-section by 6 
inches in length and provides an approximate run time of 
5 minutes at a Mach number of 2.80 and a stagnation pressure 
of 50 psia. The lower (tunnel floor) aluminum block was 
replaced by an accurate phenolic replica of the original 
lower block in order to facilitate mounting the model 
with its required tunnel floor holes while preserving the 
geometry of the original block intact. 

C. THE EXPERIMENTAL SET-UP 

Figure 4 is a schematic representation of the equipment 
arrangement. The calibrating apparatus is not shown. 

Injectant gas, manifolded from banks of four commercial 
high pressure cylinders , was supplied through a throttle 
valve into a sharp-edged orifice flowmeter. The flowmeter 
upstream pressure was measured in psia by a Heise gauge 
(Serial No. H 2946) . Differential pressure across the 
orifice was sensed by a Statham differential pressure 
transducer (Serial No. 11404) calibrated in inches of 
mercury and capable of measuring a ± 10 psi differential 
at a working pressure of 500 psig maximum. The stagnation 
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temperature of the flow was measured at the flowmeter by 
means of a chromel-alumel thermocouple. The injectant gas 
then proceeded to an expansion chamber mounted adjacent to 
the lower wall of the wind tunnel in which the flow was 
effectively stagnated. At this station the total pressure 
of the secondary flow was measured by a Heise pressure 
gauge. The total temperature of the injectant measured 
at the expansion chamber agreed closely with that measured 
downstream of the flowmeter (within 5° F) . The stagnation 
temperature at the flowmeter was taken to be the injection 
total temperature. The total temperature and total 
pressure of the primary flow was measured in the wind tunnel 
plenum chamber. 

For side force determination, the static pressures were 
measured at the pressure taps on the flat plate model by 
a Scanivalve (Model 48J4 - 1065 Scanivalve with a 0-10 psia 
transducer) calibrated in inches of mercury vacuum. The 
measurements were recorded on a Honeywell 2106 twelve 
channel Visicorder as was the pressure differential across 
the flowmeter orifice mentioned above. 

Schlieren photographs on Polaroid 4x5 type 55 P/N 
positive-negative film were taken of each run using an off- 
set Schlieren system with collimating mirrors. 

Photographs of the experimental equipment are provided 
in Figures 5 through 9. 



17 



III. THEORETICAL ANALYSIS 



From the many phenomena associated with secondary 
injection, it was chosen in this paper to deal with the 
characteristic radius which may be computed from the blast 
wave analogy model, the penetration height of the injected 
gas, the shape of the interaction or bow shock and the 
resultant normal side force developed on a flat plate. 

A. PENETRATION HEIGHT 

1 . The Blunt Body Model of Zukoski and Spaid 

The similarity of the shock shapes produced by 
secondary injection into supersonic main flows and the shock 
patterns produced by blunt axisymmetric bodies in supersonic 
flow prompted Zukoski and Spaid [Ref. 1] to develop a blunt 
body analogy to the secondary injection phenomena. Assuming 
the proper choice of a blunt body, the characteristic 
dimension of the nose (a length analogous to the penetration 
height of the secondary injectant) may be computed by equat- 
ing the drag of the required blunt body to the change in 
momentum flux of the secondary injectant. These investi- 
gators assumed: 

a. Sonic injection into a uniform supersonic flow 
with no boundary layer 

b. No mixing between primary and secondary flows 

c. The blunt body to be a quarter sphere with a 
semi-cylindrical afterbody 
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d. The injectant gases remain within the confines 
described in item C 

By assuming Newtonian flow about a body in the form of a 
quarter-sphere with a half-cylinder afterbody, Zukoski and 
Spaid obtained the force balance 

(In jection) ~ (Blunt Body) 



or 



m V.^ = (Tr/8) Y c*h^ P 
s ]X ' / ' “ P p ° 



(2) 



is the streamwise velocity of the injectant and is 
determined by assuming isentropic expansion of the injectant 
from the secondary stagnation pressure to free-stream 
pressure. C*^ is the Newtonian pressure coefficient behind 
a normal shock and is evaluated at the nose of the analogous 
blunt body. 

Solving Equation (1) for the radius of the blunt 
body (the effective penetration height h) and applying a 
nozzle coefficient of discharge correction to permit the 
use of physical nozzle dimensions in these calculations, 
the following expression was obtained: 



1 2P . Y 

h/a/H = I (ctIy?- 

°° p pt'p 



) f(74r)< 



2 (Y3-^1)/<Y3-1) 



Y +1 

^ s 



[l-(p2i) 

St 



P.,4. (Y^-1)/ (Yc.+1) 



]> 



(3) 
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2 . 



The Blunt Body Model of Cassel, Davis and Engh 



In a manner somewhat similar to the method of Zukceki 
and Spaid, work by Cassel, Davis and Engh [Ref. 2] employed 
an equivalent blunt body analogy in the development of an 
expression of effective penetration height of a secondary 
injectant. The obstacle envisioned in their model is a 
cylinder which is tangent to a flat plate, capped on the 
upstream end by a half-sphere and oriented in the streamwise 
direction of the primary flow. Assuming that the momentum 
of the secondary gases is parallel to the plate at the 
sonic line in the plume of the secondary jet, these 
investigators also equate the change in the streamwise 
momentum of the injectant with the drag of the equivalent 
body. The following equation results; 

C A*q^ - (P*-P ) A* = m [a*+V sine(90°-e)] (4) 

where A* is the projection of the sonic area on a plane 
normal to the wall, P* is the pressure at sonic conditions, 
A^ is secondary throat area, a* the speed of sound at sonic 
conditons and is the secondary nozzle exit velocity. 
Continuity implies 



m = p*A*a* = 






, 2 



for a perfect secondary gas at constant total temperature. 
Therefore, assuming that the frontal area of the injectant 
stream is effectively circular. 
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( 5 ) 



P^.* * A ^4- 2 



st 



A* 



Substituting into Equation (5) the two previous 
equations yields 

A* h 2 2 ^st . \ 

f- = a +P )^Ys[l+^sin(90°-e)]+l} 

^t ^t ^s ® ^t 

In Ref. 3, the vacuum thrust coefficient of a nozzle is de- 
fined as 



2 Yg/(Y -1) 

^f = 

o ' s 



( 6 ) 



combining the above equations, Cassel, Davis and Engh obtain 



F . 



A* = 



Y. 



V 



‘^l C^+{2/y^uJ) Yg+1 



) {r^[l+ — sin(90<’-e) ] + (-^) } 



Y +1' 
' s 



where the last two terms are approximately one for normal 
sonic injection. This reduces the expression for sonic 
injection to 



A* = F. /q, 
2o' ^1 



or 



h A* ^ 

■^ = (9-) = ( 



F. 

3° ^ 



‘3l\ 



(7) 



Crist, Sherman and Glass [Ref. 4] obtained a good fit between 
analytical and experimental results with the expression 
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or, from Equation (5) 



h^/h* = 0.71 (P*/Pj^) 



which led to the expression for e = 90' 






( 8 ) 



3 . Penetration Height as a Function of Secondary 
Mach Number 

Schetz and Billig [Ref. 5] propose a theoretical 

relationship between the values of h for varying M which 

s 

is based on the ratio of dynamic pressure at initial injec- 
tion conditions to that at mainstream conditions. For 
Y =1.4 the following relationship between penetration 
height at M =1.0 and that for M greater than 1.0 was 

5 S 



obtained 



1.2M ^ 
s 

1+0. 2 M ' 



where h^^ is the penetration height at sonic injection. 



B. SHOCK SHAPE AND SIDE FORCE PREDICTION 
1 . Blast Wave Analogy 

Hayes [Ref. 6] and Lees and Kubota [Ref. 7] proposed 
a similarity between hypersonic flow about an unyawed, 
axisymmetric body and the detonation of a line charge. In 
general, this similarity is based upon correlating the time 
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dependent pressure, density and velocity fields which are 
generated by the explosion of a line charge with the 
similar field, generated as a function of axial distance 
over the velocity of the undisturbed flow, about a body in 
a steady hypersonic flow. In their blast wave analogy 
approaches to the solution of the secondary injection prob- 
lem, Broadwell [Ref. 8] and Dahm [Ref. 9] have employed 
this similitude concept and Sakuri ' s [Refs. 10 and 11] 
solutions to the flow fields behind blast waves resulting 
from the explosion of line charges of uniform energy per 
unit length. Two basic assumptions are common to the 
development of both Broadwell and Dahm. The first is that 
the injectant gas profile is equivalent to a blunt nosed 
body and the second that the pressure field between the 
shock and a blunt nosed body is correctly predicted by 
the blast wave theory. The first assumption has been 
justified by the observed similarity of the pressure 
fields about blunt bodies to those generated by flat 
plate injection. The second is less tenable in that it 
does not hold strictly true where hypersonic flow assump- 
tions do not apply, i.e., near the nose of the blunt body 
or equivalent and at mainstream Mach numbers not in the 
hypersonic range. Both Dahm and Broadwell assume also that 
at some point not far down stream the injectant gas has 
been accelerated to the mainstream velocity. 

The energy per unit length of line charge in the 
blast wave theory may be considered analogous to the energy 
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added to the flow by the injectant. Broadwell's develop- 
ment considers this energy term to be composed of the 
secondary momentum change resulting from acceleration of 
the injectant to the mainstream velocity plus the energy 
generated by the volume addition associated with injection. 
Dahm derived this analogous energy term directly from 
thermodynamic considerations with the assumption of negli- 
gible mixing of primary and secondary streams. 

In the blast wave theory, a full 360° cylindrical 
shock wave is considered while in the case of secondary 
injection, where the flat plate is a plane of symmetry, 
only a half-cylinder shock wave is formed. Therefore, the 
analogous energy term becomes 

E = 2E' 

where E is the energy per unit length required for blast 
wave formation. 

Sakuri relates a characteristic radius , R , to the 

o 

blast wave energy and primary static pressure by: 

R = (E/27TP (10) 

o °° 

The energy term as derived by both Broadwell and Dahm is of 
the general form 

E = 2m V (11) 

s p 2 

Broadwell, by use of previously mentioned assumptions, 
ascribed to 1 O 2 the expression. 
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00 „ = 1 + ( 



l+-tL— M m T ^ 

2 ^ 

(Yp-l)M^ s pt 



( 12 ) 



whereas, from thermodynamic considerations, Dahm developed 
the form 



"2 ' 2 



Y (Y -DM ^ 

'p 'p 00 



Y 1 + — — M ^ m T . 
s , 2 «> X , p st 



(- 



.,2 ^ ^m T . ^ 

s pt 



YmT.JgM Y“1 

' o o -I- o • T-> 



Y,-l 



+ cose (13) 



Y m T . M 
'p s pt p 



2 «> 



where e is the angle between the primary flow direction and 
the injection axis. 

For normal injection, with which Broadwell was 
exclusively concerned, the expressions in Equations (12) a 
and (13) have been noted to be quite similar in form 
although their origins differ considerably in concept. 
Henceforth, 0)2 will be used here exclusively as expressed 
by Dahm in Equation (13). Therefore, 

= E/2m V 
2 ^ s p 



and Equation (10) becomes 



R = (w^m V /ttP ) 
o 2 s p' <» 



Taking 



(14) 
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R = R/R^ 

? = X/R^ 

t = x/V 

P 

where R is the radius of the interaction shock or bow shock 
at X. The shock radius predicted by blast wave theory is 
given by Sakuri as 



= Jo [1 + (^)% I t ...] (15) 

j R D D 

where J is the first order blast wave constant. A, are 
o ' k 

higher order blast wave constants, is the shock velocity, 
and a^ is sonic velocity in the undisturbed mainstream. 
Solution of Equation (15) yields 



^ J* 

R = [2aR /J ] ^ 



or, for the first order theory, 
R= [2/J^'^] " (x/M^)'^ 



and, for the second order theory. 






2J 

o o 

The theoretical bow shock shape for any 
tion parameters can thus be determined. 

In Broadwell's development, the 
is described by 



(16) 

given set of injec- 
interaction force 
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F = 2 (P-P ) drdx 

w o o 



where 



P = P [P^,J^,M^,x,g(n)] 
which yields 



F = F [E,J ,L,g(n)] 
w ° 



where g(n) is a non-dimensional pressure ratio. 

In a similar fashion, Dahm developed the expression 
for second order theory interaction side force 



12u, 



3/2 



P„A^ 



= (- 



TT^J 



M 2 % 

3/2 ' M ^ 



m V 3/4 

-S-E) 

P^.A^'^ 



o 

This expression applies only at small values of the blast 
wave analogy parameter s/R^M^ (less than 0.0654 according 
to Cassel, Davis and Engh) . For larger values of this 
parameter, i.e., for lower mass flow rates, the following 
relationship was developed in Ref. 9 for circular port 
injection 



F. = (2/iT)co^m V M z 
1 ^ 2 s p <» 



where 

z 





J k., 
o 2 



k^b 



3/2 






J k 
o . 



k^b 



3/2 



) in ( 



4- + y l+t 

b>5 / b 



) J 



(17) 
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in which b = 4> = t>R when R is evaluated at s. 

Reference 9 evaluates = 0.4973 and = -1.103 for 

Yp = 1.4. 

2 . Upstream Model 

Consideration of the secondary injectant as being 
analogous to a blunt nosed obstacle in the primary stream 
lead to the two-dimensional upstream model of Vinson, 

Amick and Liepman [Ref. 12]. Here the weak shock formed 
by boundary layer separation in the high pressure region 
ahead of the secondary stream or the analogous blunt body 
yields the equivalent penetration height of the secondary 
injectant. The penetration height thus determined, the 
boundary layer detachment angle, the mainstream pressure 
and the pressure behind the shock in turn yield the magni- 
tude of the interaction side force developed. 

Wu, Chapkis and Mager [Ref. 13] extended the pre- 
vious model to three dimensions by assuming that the 
separation region ahead of the secondary nozzle is conical. 
Their model of the interaction phenomena leads to the 
expression : 

2 

F./F . = {cot 6 + tan (a+e' ) } /2 y /2/(v +1) 

1 3 TT S ' S 

h 

.{(^ - l)/[2Yg + (Yg-l) (^ - 1)]} (18) 

where a is the primary nozzle half angle and is zero for the 
purposes of this paper, e' = 90-e, P is the pressure in the 
area of the detached region ahead of the bow shock, 6 the 
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boundary layer attachment angle and 6 the associated conical 
shock angle. 

3 . Downstream Model 

Walker, Stone and Shandor [Ref. 14] propose a model 
in which the immediate mixing with the mainstream gas of 
a low flow rate injectant in a constant area is suggested. 
After mixing, the gas is assumed to expand isentropically 
to the mainstream pressure thus inducing a side force. 

The following expression for side force is obtained: 



F. 

1 





[2(M^^-1)]'^ 








Y. 



Y +1 
s 



) ( 



m T .h 

s pt 

m T . 
p st 



) G 



(19) 



where 
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IV. DESCRIPTION OF THE EXPERIMENT 



A. PARAMETERS 

The secondary nozzles constructed for use in this 
experiment were of Mach numbers 1.00, 1.73 and 2.13. A 
description of the calibration of the supersonic nozzles 
is contained in Appendix A. The throat diameter of all 
nozzles was a nominal 0.100 ± 0.001 inches. Nominal 
inclination angles of 90°, 80°, 70°, 60° and 50° were 
tested at each injection Mach number. 

During all runs, the primary flow Mach number remained 
constant at M^ = 2.80. The primary total pressure was a 
constant 50.0 psia. In successive closely spaced runs, 
the primary total temperature decreased slightly as the 
tunnel walls were cooled by the flow. This temperature 
decrease was slight, having an overall variance of ± 15° F 
from the mean total primary temperature for all runs. For 
purposes of the interpretation of the results of the 
theoretical computations involving the primary flow 
parameters, this temperature variation was considered 
negligible . 

The secondary total pressure was increased by increments 
of 25.0 psia from 25.0 to 250.0 psia providing a series of 
ten secondary-to-primary total pressure ratios for each of 
the 15 nozzles tested. The secondary gas reservoir total 
temperatures varied only slightly due to minor fluctuations 
in the ambient temperature over the period of the testing. 
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Therefore, the total temperature of the secondary flow after 
it was throttled to the desired total pressure also varied 
slightly. This effect was also neglected for purposes of 
interpretation . 

B. EXPERIMENTAL PROCEDURE 

The test runs were conducted in ascending order of 

and at a given secondary Mach number in descending order 

of e. Ten runs were conducted at each nozzle configuration 

at the previously described total pressure ratios using 

nitrogen as the secondary gas. Two nozzle configurations 

(M = 1.0, e = 90° and M = 1.73, e = 68°) were tested 
s s 

similarly with argon and helium as the secondary injectant. 

Individual runs averaged less than one minute in 
duration. During each run the following observations were 
made : 

1. The static pressure at each of the 47 flat plate 
pressure taps v/as measured by a Scanivalve and 
recorded on a Visicorder. 

2. The differential pressure through the flowmeter 
was measured by a differential pressure trans- 
ducer and recorded on a Visicorder. 

3. The upstream flowmeter pressure was measured by 
pressure gauge. 

4. The injection total pressure was measured by 
pressure gauge. 
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5 . 



The primary total pressure was measured by pressure 
gauge . 

6. The total temperature of the secondary flow 
through the flowmeter was sensed by a chromel- 
alumel thermocouple and displayed on a millivolt- 
meter . 

7. Primary total temperature was sensed by a chromel- 
alumel thermocouple and displayed on a millivolt- 
meter . 

8. A Schlieren photograph was taken of the test 
section . 

As described earlier, the total temperature of the secondary 
flow at the flowmeter was taken to be the injection total 
temperature . 

Tare pressures were observed at the plate surface with 
the tunnel operating but without secondary injection. 

These were then subtracted from the respective pressures 
measured during injection to determine the interaction 
pressure rise. This pressure rise was then integrated 
over the plate surface for determination of the interaction 
side force. Tare observations were made prior to each of 
the series of 10 runs for a given nozzle configuration. 

In all cases, tare pressures were observed with a nozzle 
for normal injection installed. 

Penetration heights and bow shock coordinates were 
measured by projecting the negatives of Schlieren photo- 
graphs onto a measuring grid. The grid was constructed to 
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twice the actual scale of the model and reference points 
were provided to insure that the grid was properly 
oriented with respect to the projection of the negative. 

In order to determine the effect of the protrusion of 
the downstream side of the nozzles at e < 90°, a duplicate 
sonic normal nozzle was constructed. A dummy obstruction 
representing the protrusion of an inclined nozzle was then 
inserted just downstream of the normal sonic nozzle orifice. 
This nozzle was tested in order to determine the effect 
of the obstruction on the pressure distribution over the 
plate and on penetration height. 

A normal sonic nozzle was constructed with a cylin- 
drical constant area throat section 0.01 inches in length 
and tested for comparison with the sharp-edged-throat 
nozzles used in this experiment 
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V. EXPERIMENTAL RESULTS 



A. PENETRATION HEIGHT 

Zukoski and Spaid have stated that the height of 
penetration of the secondary gases into the main flow is 
a dimension characteristic of the scale of the interaction 
disturbance. This parameter, non-dimensionalized by the 
nozzle throat diameter and the square root of the nozzle 
discharge coefficient (i.e., by the effective nozzle 
throat diameter) has been utilized by previous investiga- 
tors [Refs. 1, 15 and 16] to correlate injection penetration 
height with secondary-to-primary total pressure ratios for 
normal, sonic injection. 

In this experiment, the discharge coefficient was found 
by dividing the actual flow rates through the nozzles by 
the corresponding mass flow rate calculated for ideal gas 
assuming isentropic flow. The actual flow was measured by 
means of a sharp-edged orifice flowmeter placed in the injec- 
tant supply line. The flowmeter was constructed and the 
flow through it evaluated in accordance with Ref. 17. 
Discharge coefficients thus obtained are shown versus the 
inverse square root of nozzle throat Reynolds number in 
Figures 10, 11, and 12. These results are in general 
agreement with the discharge coefficients of similar nozzles 
determined by Amick and Hays [Ref. 16] and by Chrans and 
Collins [Ref. 15]. 
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As described in detail earlier, measurements of the 
penetration heights were taken from the Schlieren photo- 
graphs which were made of each run. This dimension was 
measured along a normal from the surface of the flat plate 
to the point of highest penetration of the secondary flow. 

For injection at a given Mach number, penetration height 
correlations in non-dimensional plots of (h/d/c) sin e versus 
(Pg^/Pp^) sin e for all e were obtained. The results are 
shown in Figures 13, 14 and 15. For the special case of 
normal (sin e = 1) , sonic injection Figure 13 is equivalent 
to a plot of h/d/c as a function of P ./P . which, as 

St pt 

stated earlier, has seen previous use. The present 
application extends the generality of the earlier correlation. 

From Figures 13, 14 and 15 it can be seen that 

p 

( ) sin e = [ (^ — ) sin e]^ 

d/c pt 

where the exponent 3 varies with Mach number. In Figure 
13, 3 was found to be approximately equal to 0.53. This 
was in close agreement with the value of 3 = 0.5 in the 
theoretical development of Zukoski and Spaid for normal 
sonic injection. As can also be seen in Figure 13, the 
observed data agrees closely with the theoretical pene- 
tration height for normal, sonic injection of Cassel, Davis 
and Engh. The Riemann disc height as a function of pressure 
ratio for an underexpanded sonic nozzle exhausting into 
still air as determined in Ref. 18 is also shown in Figure 



35 



13. This curve agrees closely with both the observed and 

the theoretical correlations. 

As expressed by Schetz and Billig in Equation (9), the 

penetration height is expected to increase as a function 

of Mach number. In Figures 13, 14 and 15 this is borne out 

in a qualitative sense. Figure 16 is a reproduction of the 

mean lines defined by the data in the preceding three 

figures. As shown in this figure, the penetration height 

does increase with increasing secondary Mach number. 

However, values of the exponent 3 are also seen to increase. 

As indicated by the theoretical plots also shown in Figure 

16, this was not predicted by Schetz and Billig. At 

M =1.73 the observed value of 3 is approximately 0.62 and 
s 

at M =2.13 approximately 0.63. 
s 

B . SHOCK SHAPE 

Dahm and Broadwell suggest that the characteristic radius 
is a dimension representative of the scale of the inter- 
action. Dahm has shown theoretically that a non-dimensional 
plot of x/R^ = X versus = R provides correlation of 

bow-shock coordinates. Experimental verification of this 
correlation for normal sonic injection was obtained by 
Chrans and Collins with variations in mass flow rate and 
secondary-to-primary total temperature ratio. The present 
study has shown that these coordinates are not sufficiently 
general. Simultaneous correlation is not achieved when data 
is collected at various values of inclination angle e. As 



36 



shown for nitrogen injection in Figures 17 through 30, it 
was, however, found that bow shock shape correlated in a 
plot of X versus R for any specific e. Correlation with 
respect to all values of epsilon was obtained by plotting 
shock shapes in the non-dimensional coordinates x sin e 
versus R sin e. For normal injection, where sin e = 1, the 
coordinates in effect revert to x and R. Bow shock shapes 
are correlated in this manner in Figures 31, 32, and 33 for 
nitrogen, argon and helium injected at a secondary-to- 
primary total pressure ratio of 3.0. It can be seen from 
the previous data correlations for various pressure ratios 
in Figures 17 through 30 that the results here are repre- 
sentative of all pressure ratios. 

The data scatter in Figures 31, 32, and 33, although 
coherent, is separated to a slight degree with respect to 
secondary Mach number. In general, the data representing 

injection at M =2.13 lie along the upper bound of the 
s 

scatter band, that for M = 1.73 in the center and the data 

o 

for sonic injection along the lower bound. The parameter 

R , by which bow-shock coordinates were non-dimensionalized 
o 

in these figures, was calculated from Equation (14) in 

which R^ is a function of a >2 • As can be seen from Equation 

(13), ( 0 - is a function of M only for e < 90°. This depen- 

dence on is sufficient to cause shock shapes to correlate 

very closely for all e at a given M . Correlation at 

s 

differing values of M is, however, not so precisely 

s 

obtained. In the previous section, h was suggested as a 
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characteristic dimension representative of the scale 
of the disturbance. In Refs. 1 and 15 bow shock shape 
correlation for the case of normal, sonic injection was 
obtained by plotting x/h = x' versus R/h = R' . Corre- 
lation of shock shapes for all inclination angles was 
obtained in Figure 34 by plotting R' sin e versus x' sin e. 
For normal injection (sin e = 1) these coordinates are 
identical with the coordinates previously used for correla- 
tion of shock shapes formed by normal sonic injection. It 
was found that when correlated in this manner, the data 
for various secondary Mach numbers was randomly distributed 
throughout the scatter band rather than forming obvious sub- 
groups by Mg within it as was the case when shock coor- 
dinates were non-dimensionalized by R . 

The mean lines representative of the shock shape cor- 
relation data in Figures 31, 32, and 33 (for nitrogen, 
argon and helium) are reproduced for comparison in Figure 
35. The dependence of shock shape on molecular weight 
previously noted by Chrans and Collins is seen. At a 
given value of x sin e the non-dimensional radius, R sin e 
increased with increasing secondary molecular weight. Blast 
wave theory does not indicate such a dependence. It has 
been suggested [Ref. 15] that the mixing of the secondary 
gases with the primary stream (assumed to be negligible in 
the development of the blast wave analogy) may account for 
this dependence. 
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When plotted in terms of x' sin £ versus R' sin e, the 



same molecular weight dependence is shown. 



C. SIDE FORCE 



An experimental examination of the interaction side 
force developed by secondary injection was also conducted. 
Side forces on the flat plate were calculated from the 
interaction pressure rise on the plate surface as detailed 
earlier. The interaction side force was calculated for all 
runs . 

In this study, side forces are presented in the form 
of a non-dimensional force amplification factor A. The 
amplification factor is defined as 



A = 



F . + F . sin e 



F . son e 
I 



which relates the interaction side force to the component 
of secondary nozzle thrust normal to the plate surface. 
Reproducibility of the observed values of A were consis- 
tently good. Figure 36 (a) shows representative data 

scatter in a plot of A versus primary-to-secondary total 
pressure ratio. 

In order to prevent disruption of the symmetry of the 
secondary flow by the physical hardware of the model, all 
nozzles were constructed with the upstream edge of the 
nozzle exit planes coincident with the plate surface. As 
a result, the downstream edge of the nozzles protruded 
slightly into the primary flow. It was necessary to 
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determine the effect of this protrusion. A sonic, normal 
nozzle was constructed with an obstruction representing an 
inclined nozzle installed immediately downstream of the 
nozzle orifice. This nozzle was then tested for comparison 
with the normal sonic nozzle used in the experiment. The 
values of A at various injection pressure ratios thus 
obtained agreed within the normal experimental scatter 
encountered. A normal sonic nozzle having a short (0.01 
inches) cylindrical, constant area throat (square edged) 
was also tested for comparison with the sharp edged nozzles 
used in this experiment. The amplification factors were 
found to agree within the bounds of anticipated scatter. 

The results of these comparisons are shown in Figure 36 (b) . 

Schlieren photographs showed no difference in either bow 
shock shape or penetration height in these comparisons of 
nozzle configuration. 

In the development of their blunt body models, Zukoski 
and Spaid, and others [Refs. 2 and 8] suggest that the 
characteristic length is a function of the change in 
momentum experienced by the secondary gases. 

As has been shown earlier, penetration height increases 
as epsilon decreases (i.e., as the upstream component of 
momentum increases) when all other parameters are held con- 
stant. It is, therefore, to be expected that as the scale 
of the disturbance increases the interaction side force 
will increase accordingly. In Figure 37, at a repre- 
sentative secondary-to-primary total pressure ratio of 
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3.0, a plot of interaction side force versus epsilon is 

shown. It may be seen that a decrease in epsilon is 

accompanied by an increase in observed side force. This 

figure also shows a decreasing trend in interaction side 

force as M increases. The interaction side force for 
s 

sonic injection and that for injection at = 1.73 at 

various values of e are nearly equal. For injection at 

M =2.13, however, a decrease in observed side force 
s 

occurs. This decrease is not predicted by Dahm. Reference 
5 suggests that the gains which theoretically appear to be 
associated with supersonic injection may be countered by 
the losses involved in turning a supersonic secondary jet 
downstream. The turning of a sonic secondary jet is 
largely shock free and the losses are thus small compared 
to the turning of a supersonic jet. 

In Figures 38, 39, and 40, interaction side force is 
plotted in terms of the nondimensional force amplification 
factor A versus secondary-to-primary total pressure ratio. 
The amplification factor A, at a given value of pressure 
ratio, is seen to increase with decreasing epsilon. Part 
of this increase is due to the decreasing value of the 
normal component of jet thrust (F^ sin e) as e decreases 
and part is due to the increasing value of F^ as e decreases 
which is shown in Figure 37. 

It is noted in Ref. 19 that an empirical correction 
factor of 0.51 when applied to Dahm's second order theory 
values of side force provides a semi-empirical value of 
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side force very nearly comparable with previous experimental 
data. The shaded areas on Figures 38, 39 and 40 represent 
the region of the value of A computed from this semi- 
empirical correction to Dahm's theory. It can be seen 
that considerable data scatter exists in the area of low 
injection pressure ratios and that observed and semi- 
empirical values differ considerably. Conversely, at the 
higher pressure ratios, scatter decreases and agreement 
with the semi-empirical values of A improves considerably. 
Also, as Mg increases, a comparison of Figures 38, 39 and 
40 shows that overall data scatter decreases and overall 
agreement with the semi-empirical values of A increases. 

In Figures 41, 42 and 43, the data shown in Figures 

38, 39 and 40 is presented in a more general manner. The 

previously observed successes in correlation of parameters 

suggest, as e is varied, that all features of the flow 

scale by sin e when properly non-dimensionalized. A degree 

of correlation was achieved in Figures 38,39, and 40 in 

a plot of A sin e versus nigVj sin e. Again, at lower values 

of nigVj sin e(i.e., lower injection pressure ratios) data 

scatter is pronounced and agreement with empirically 

corrected theory is not good. As m V . sin e increases, 

s 3 

however, data scatter decreases and agreement with the 
corrected theory continually improves. As increases, 
comparison of Figures 41, 42 and 43 also shows considerable 
improvement in the degree of correspondence of data with 
the curve representing the semi-empirical values of A. 



42 



The degree to which observed values of side force agree 
with the side force actually developed is a function of the 
fineness of the grid spacing of the pressure taps. It 
was found in this experiment that both data scatter and 
disagreement with corrected theory were pronounced at the 
lower values of injection pressure ratio (or the related 
lower values of sin e) . In this area of operation, 

the characteristic length h was of the order of or less 
than the average distance between pressure taps. As h 
increased to a value of the order of twice the pressure 
tap spacing or greater (i.e., for values of (Pg^/Pp^.) sine 
greater than 2.0) data scatter decreased and agreement with 
semi-empirically predicted values of A increased. Above 
values of (P^^/Pp^) sin e > 2.0 or the corresponding values 
of sin e > 25, the data in Figures 38 through 43 is 

considered reliable. In this area, the data in Figures 
41, 42, and 43 shows a consistant degree of correlation 
which agrees well with expected values. 

As determined by Zukoski and Spaid, the surface area 
ahead of the injector is subjected to the highest inter- 
action pressures developed during secondary injection. 

It was observed in this experiment that the pressures 
developed downstream of the injector were generally of 
little significance in calculating the value of inter- 
action side force. The model of the interaction proposed 
by Wu, Chapkis and Mager depends exclusively upon flow 
phenomena upstream of the nozzle for determination of the 
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interaction side force. A limited investigation of this 
theory was undertaken. The results of the comparison of 
observed data and that calculated according to Equation 
18 is shown in Figure 44 in a plot a A versus 
The calculated values were determined by measuring the 
boundary layer detachment angle 6 from the Schlieren photo- 
graphs of the respective runs. These measurements were 
accomplished on an optical comparitor. The shock angle 0 
was not shown on the photographs with sufficient clarity 
to permit measurement. Therefore, the value of 6 corres- 
ponding to a cone of half-angle 6 was taken from Ref. 20. 
The pressure behind the conical shock was also determined 
from Ref. 20. The results of the computations thus 
achieved are shown by the curve in Figure 44. The broken 
portion of the curve represents areas in which photographic 
quality did not permit accurate measurement of 6. Observed 
and calculated values agreed closely. 

The side force determined by the downstream theory of 
Walker, Stone and Shandor was evaluated in accordance with 
Equation 19 for normal sonic injection. The value of A 
at low values of pressure ratio was found to be of the 
order of 6 , in considerable disagreement with observed and 
previously computed values. 

In Figures 45 and 46, values of A are plotted against 

pressure ratio and mV. sin e, respectively, for nitrogen, 

s D 

argon and helium for two different nozzle configurations. 
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within the previously encountered range of data scatter, 
no dependence of A on the molecular weight of the injectant 
was observed. 
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VI. CONCLUSIONS 



Penetration height increases both with an increase in 
and with a decrease in e , all other parameters being 
held constant. Correlation of penetration height for a 
given value of is obtained in a plot of (h/d/c) sin e 
versus normal sonic injection these 

coordinates become the familiar h/d/c versus 
the case of sonic injection , the slope of the correlation 
curve is approximately 0.5, i.e., 

h ^4- 

sine®[ (p^) sin e ] 

d/c pt 

or for the case of normal sonic injection 

0.5 ^ 0.5 

h/d/c 

This is in close agreement with the theories discussed. 

As increases, the value of the exponent increases. At 

M =1.0 the value of the exponent is approximately 0.53, 
s 

at M =1.73 approximately 0.62 and at M =2.13 approx- 
s s 

imately 0.63. 

Non-dimensional bow shock shape plotted in the coordinate 
set (x, R) remains essentially constant with varying . 

The non-dimensional radius at a given value of x (where x 
is measured downstream from the nozzle centerline) increases 
with decreasing e. Bow shock shapes correlate in non- 
dimensional plots of X sin e versus R sin e or x' sin e 
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versus R' sin e. The correlation curves are coincident 
with those for normal sonic injection plotted in the coor- 
dinates X versus R or x' versus R' . 

The empirical correction factor of 0.51 applied by 
Dahm to the values of interaction side force calculated 
from his second order theory was qualitatively substan- 
tiated. Correlation of side force data in a plot of A 

sin e versus m V. sin e is indicated for alls and M . 

s j s 

In order to obtain valid side force data by means of inte- 
grating interaction pressure rise over a flat plate, the 
spacing of the pressure taps should be on the order of or 
less than one-half the smallest penetration height expected. 

Dahm's theory predicts an increase in the magnitude of 
side force with decreasing e, all other parameters remain- 
ing constant. This was qualitatively confirmed. Side 
force at a given value of e was observed to decrease with 

increasing M . Dahm's theory predicts no change with 
s 

increasing for normal injection and a slight increase 
in side force with increasing as e decreases. 

For the same values of injectant momentum, the 
experimental results indicate, within the degree of data 
scatter expected, that side force does not vary with changes 
in the molecular weight of the secondary gas. At the nozzle 
throat, i.e., at choked flow, the secondary mass flow rate 
is reduced as the molecular weight of the injectant de- 
creases. It is beneficial, therefore, to inject gases of 
low molecular weight in order to minimize the total weight 
of secondary fluid which must be carried by the vehicle. 
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APPENDIX A 



SUPERSONIC NOZZLE CALIBRATION 



The supersonic nozzles utilized in this experiment were 
constructed with 10° divergence half-ang]es on either side of 
a sharp edged throat. The diameters of the throats of all 
nozzles were 0.10 ± .001 inches. Throat diameters were held 
to this tolerance in the following manner. Two carefully 
manufactured rods were constructed, one having a diameter 
of 0.10 - .001 = 0.099 inches and the other having a diaiteter 
of 0.10 + .001 = 0.101 inches. The 10° half-angle approach 
portion of the nozzle was drilled out. The portion of the 
nozzle downstream of the throat was then drilled (also with 
the 10° half-angle bore) until the smaller diameter rod but 
not the larger diameter rod was free to pass through the 
nozzle throat. Thus, the throat diameter was within the 
diameter difference of the testing rods, or 0.10 ± .001 
inches . 

The nozzle having the lowest supersonic Mach number was 

constructed to a nominal area ratio of 1.6900 corresponding 

to M = 1.998 in Ref. 20. The diameter tolerance allowed 
s 

at the throat was also incurred at the exit plane of the 
nozzle. The area ratios of the nozzles constructed, there- 
fore, varied between 1.6823 and 1.6979. This causes a Mach 

number variation between M = 1.996 and M = 2.001 for one 

s s 

dimensional, inviscid isentropic flow calculations,. 
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The higher Mach number nozzle was constructed to a 

nominal area ratio of 2.6569 corresponding as above to 

= 2.508. Tolerances in manufacture caused a possible 

spread in area ratio between 2.6366 and 2.6777 or, in 

Mach number, between M = 2.500 and 2.516. 

s 

The nozzles, however, are actually three dimensional. 
Additionally, the viscosity of the secondary gas will 
cause boundary layer effects to be present, especially 
downstream of the nozzle throat where a turbulent boundary 
layer will exist. The nozzle in effect will sense a 
smaller area ratio than that for which the nozzles were 
designed. The overall effect is to decrease Mach number 
observed from that determined by one dimensional, isen- 
tropic, inviscid calculations. 

Calibration was accomplished by suspending a thin 2.5° 
half-angle cone (with its apex pointed upstream) in the 
flow exiting the nozzles. The apex of the cone was 
centered in the exit plane of the nozzle. Schlieren 
photographs were then taken of the flow over a range of 
nozzle flow stagnation-pressure-to-atmospheric-pressure 
ratios . 

The conical shock angles, 0, were measured from the 
photographs on an optical comparitor. The angles thus 
determined were 6= 35.5° ± 0.5° for the low Mach number 
supersonic nozzle and 0 = 28.0° ± 0.5° for the higher Mach 
number supersonic nozzle. These correspond, respectively. 
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2.13 ± 0.05 from the conical 



to M = 1.73 ± 0.05 and M = 
s s 

shock curves of Ref. 20. 

Taking the actual uncertainty in caused by variation 
of the physical dimensions in the nozzles to be of the same 
magnitude as that determined by one-dimensional flow 
calculations, the nozzles constructed operated at = 1.73 
± .002 and 2.13 ± .008. In comparison with the degree of 
calibration uncertainty present, the dimensional variation 
allowed in nozzle construction is, therefore, insignificant. 
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Injector 





FIGURE I 

EXPLODED VIEW OF MODEL 
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Injector Assembly 
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FIGURE 2 

TOP VIEW OF FLAT PLATE 
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FIGURE 3. TYPICAL NOZZLE DESIGN scale = 3:i 
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FIGURE 4 

EQUIPMENT SCHEMATIC 
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FIGURE 8. Experimental Equipment 
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NON-DIMENSIONAL PENETRATION HEIGHT 
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FIGURE 13 



NON-DIMENSIONAL PENETRATION HEIGHT 
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FIGURE 14 



NON-DIMENSIONAL PENETRATION HEIGHT 
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FIGURE 15 



NON-DIMENSIONAL PENETRATION HEIGHT 
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FIGURE 16 




FIGURE 17 
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FORCE AMPLIFICATION FACTOR VS. TOTAL PRESSURE RARO 

(a) ILLUSTRATES TYPICAL DATA SCATTER 
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FORCE AMPLIFICATION FACTOR (Ms =1.00) VS. NORMAL MOMENTUM FLUX 
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FORCE AMPLIFICATION FACTOR ( Ms=2.l3) VS. NORMAL MOMENTUM FLUX 
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FIGURE 44 
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FORCE AMPLIFICATION FACTOR VS. TOTAL PRESSURE RATIO 
A COMFARBON FOR GASES OF VARIOUS MOLECULAR WEIGHT 
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